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Abstract — In this paper, a robust nonlinear control scheme 
is proposed for a nonlinear multi-input multi-output (MIMO) 
system subject to bounded time varying uncertainty which 
satisfies a certain integral quadratic constraint condition. The 
scheme develops a robust feedback linarization approach which 
uses standard feedback linearization approach to linearize 
the nominal nonlinear dynamics of the uncertain nonlinear 
system and linearizes the nonlinear time varying uncertainties 
at an arbitrary point using the mean value theorem. This 
approach transforms uncertain nonlinear MIMO systems into 
an equivalent MIMO linear uncertain system model with 
unstructured uncertainty. Finally, a robust minimax linear 
quadratic Gaussian (LQG) control design is proposed for the 
linearized model. The scheme guarantees the internal stability 
of the closed loop system and provides robust performance. 
In order to illustrate the effectiveness of this approach, the 
proposed method is applied to a tracking control problem for 
an air-breathing hypersonic flight vehicle (AHFV). 

I. INTRODUCTION 

The general problem of feedback linearization for nonlin- 
ear systems with uncertainty has been typically approached 
in the literature by imposing some conditions on the un- 
certainty description which are known as matching condi- 
tions [1] and the strict triangularity condition [2]. Methods 
considering mismatched uncertainties also exist, in which 
uncertainties are decomposed into matched and mismatched 
parts. These methods typically require the mismatched parts 
not to exceed some maximum allowable bound [3]. In 
an attempt to solve a related issue arising in feedback 
linearization, in our previous work [4], [5] we propose a 
method of robust feedback linearization to feedback linearize 
a nonlinear system with uncertainties by representing the 
uncertainties in a realistic way and relaxing the matching 
condition requirements on the description of the uncertain- 
ties. In this approach, we linearized the nominal part of 
the system using the feedback linearization approach and 
linearized the remaining nonlinear terms with respect to each 
uncertainty and state to obtain an acceptable linear form for 
the uncertainty model at arbitrary operating points. However, 
in order to express the system in a more convenient set 
of coordinates, we have defined a diffeomorphism T which 
depends on the nominal values (without uncertainty) of the 
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parameters. This definition of the diffeomorphism requires 
that the system either satisfies the generalized matching 
conditions [1] which are relaxed versions of the matching 
condition or allows for additional uncertainty inputs in the 
system. 

In order to relax the generalized matching condition and 
the strict triangularity requirement, in this paper, we have 
introduced a notion of an uncertain diffeomorphism. This 
definition of uncertain diffeomorphism is similar to the 
one used in our previous work (see [6]). Furthermore, in 
order to deal with the nonlinear uncertain terms which are 
subject to time varying uncertainty, here we use a mean 
value approach similar to the approach used in [5]. The 
uncertain diffeomorphism used in this paper is the function 
of system states and uncertain time varying parameters. 
Generally, feedback linearization of higher order systems 
involves higher order derivatives of the system's outputs 
being required to be measurable. However, in the real world 
applications, it is not possible to measure or manipulate all 
of the output derivatives (new states) algebraically especially 
in the presence of uncertain diffeomorphism. The minimax 
LQG control approach solves this problem by using output 
feedback and estimating the unmeasured states in the pres- 
ence of the uncertainties. The main idea in our approach 
is to suppress the perturbations arising from the nonlinear 
uncertainties using a minimax LQG controller [7]. 

In the later part of this paper, we apply the proposed 
method to design a tracking controller for velocity and 
altitude tracking of an air-breathing hypersonic flight vehicle 
(AHFV) in the presence of input coupling and flexible mode 
effects. In this paper, we solve the AHFV tracking control 
problem by designing a robust minimax LQG controller 
in combination with the robust feedback linearized model 
proposed in this paper In the minimax LQG scheme, the 
uncertain states are estimated by using a robust Kalman filter 

The main contributions of this work as compared to 
previously published work [4], [5], [6] are as follows: 

1. Feedback linearization of uncertain systems subject to 
time varying uncertainty using an uncertain diffeo- 
morphism along with a mean value approach. 

2. Estimation of uncertain states and design of a tracking 
controller using the minimax LQG design method for 
the linearized model. 

The paper is organized as follows. Section describes 
the class of nonlinear systems and uncertainties considered 
in the paper A complete derivation of the robust feedback 
linearization of the uncertain system is presented in Sec- 
tion |III] In Section |IV] the minimax LQG control design 



method is presented for the feedback hnearized system with 
Unearized uncertainty. For the case study using an uncertain 
nonhnear model of the AHFV, the uncertainty modehng and 
control design methods with tracking simulation results are 
presented in Section |V] Conclusions are presented in Section 

m 

II. System Definition 

Here, we consider an uncertain multi-input multi-output 
nonlinear system having same number of inputs and outputs 
and which is subject to time varying uncertainty p{t): 



for the sake of brevity): 



x{t) = f(x(t),p{t)) + Y,9k{xit),pit))ukit), 

k = l 

yi{t) = i'i(a;(t)), i = 1,2, ■ • ■ ,m 



(1) 



where x{t) G R", u{t) = [ui Umf € R™ and y{t) = 

[yi----ym]'^ G R™- The nonlinear functions f{x{t),p{t)), 
and gk{x{t),p{t)) and i'i{x{t),p{t)) for i = 1, • • • ,m are 
infinitely differentiable (or differentiable to a sufficiently 
large degree) functions of their arguments. Also, p{t) G R'' G 
O is a vector of unknown parameters or disturbances which 
takes values in the set C R^. The subscript indices k 
and i indicate kth and ith elements of the corresponding 
vectors respectively. The full state vector x{t) is assumed 
to be available for measurement and the uncertainty in the 
system satisfies an integral quadratic constraint condition 
(IQC) (see [7]). It is assumed that the uncertain functions can 
be written as f{xit),p{t)) = fix{t),po) + ^f{x{t),p{t)) 
and g{x{t),p{t)) = g{x{t),po) + Ag{x{t),p{t)) where po is 
the nominal value of the parameter. In addition, the uncertain 
functions Af{x{t),p{t)) and Ag{x{t),p{t)) are smooth and 
contain all the uncertainties in the system, including dis- 
turbances and uncertain nonlinear terms. Furthermore, there 
exist an isolated equilibrium point which is not affected by 
the vector p{t); i.e. /(O) = 0, and A/(0,p(t)) ^ and 
system has full relative degree with respect to the regulated 
output. 

III. Robust Feedback Linearization 

In this section a robust feedback linearized method is used 
to linearize ([T i using a technique developed in our previous 
work; see [4], [5]. We decompose the system ^ into nominal 
and uncertain parts as follows: 

= /o(2:(t),po) + ^ 5fco(a;(t),po)ut(t) 



Nominal pan 



(2) 



Uncertain pml 

yi(t) = i^i{x{t)), 1 = 1,2, ■■■,m. 

The nominal nonlinearities in the equation ([T]l can be 
canceled using a standard feedback linearization approach 
[8]. Let us assume that for i — 1,2, ■ ■ ■ ,m is the relative 
degree of each regulated output, the Lie derivative of each 
output i^i, ri number of times, for each subsystem can be 
written as follows (we drop the argument t from the functions 



(3) 



— 1 

k=l '- 

In order to write the system in a form suitable for feedback 
linearization, we write the rf' derivative of each output as 
follows: 



= f»(x,po) + g,(x,po)u 



(4) 



where. 



g,{x,po) = 



'fo 



r ^2 - 



and the Lie derivative of the functions Vi with respect to 
the vector fields / and gk are given by 



LfVi ='^^^j^f, L^jVi = Lf{L'j ^Vi{x)), Lg^(vi) = 



dui {x) 
dx 



-gk- 



The nominal feedback linearizing control law 

u = -g,{x,poy^f*{x,po) + g^{x,po)~^v (5) 

partially linearizes the input-output map (IDl in the presence 
of uncertainties as follows: 



Vl 






+ 







AW['{x,u,po,Ap{t)) 
AW:r(x,u,po,Ap(t)) 



(6) 



Nominal part Uncertainty piwl 

where AWp {x,po, Ap{t)) = 



and 



E ^Aofc [Lf-ii^t)]uk, y* 
k=i " 

V = [vi....v„i]'^ is the new control input vector. Furthermore, 
we define an uncertainty vector AWi which represents the 
uncertainty in each derivative of the z*'* regulated output as 



AW^{x,u,po,Ap{t)) = 



AWl{x,u,po,Ap{t)) 
AW^{x,u,po,Ap{t)) 

AWp{x,u,po,Ap{t)) 



(7) 



and write yi for i = 1,2, ■ ■ ■ ,mas given below. 



r yl 1 
y^ 




- - 
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AW,\x,u,po,Ap{t)) 
AW,^x,u,po,Ap{t)) 


■ (8) 










_ AWp{x,u,po,Ap{t)) _ 





Let us define an uncertain diffeomorphism for each partially 
linearized system in ([8]i for z = 1, • • • , m as given below: 



(9) 



X^ = Ti{x,p{t)) = 

[ jvi- Vic Vt - Vic yj ■■ y[" 

Using the diffeomorphism (|9]) and system (O we obtain the 
following: 

X = Ax + Bv + AW{x,v,po,Ap{t)), (10) 

IT 



where x{t) 

[Vi_ V2-- - Vr 

AWix,v,p{t)) = [AWii-),AW2i-), ■■■ ,AW_, , 
a transformed version of AW{x,u,p{t)) and AWi{-) 

^MT f_ o - 1,2, ••• , TO. Also, 



[xi{t),--- ,Xmit)Y e M«, V{t) = 
G is the new control input vector, 



[0,0, 



A 







for i 




An 



Bi 











In order to obtain a fully linearized form for ( fTOb . here, 
we use a similar approach as used in [5]. In this work, we 
perform the linearization of AW(x,v,p) using the gener- 
alized mean value theorem [9], [10] such that no higher 
order uncertain terms exist after the linearization process. 
Since in this scheme an uncertain diffeomorphism is used, 
therefore this scheme provides a bound which would be less 
conservative than the bound obtained in [5]. 



Theorem 1: [10] Let w. 



be differentiable on 



with a Lipschitz continuous gradient Vw^ . Then for 
there is a c = x + t{x ^ x(0)) with 



given X ™d x(0) in 
i e [0, 1], such that 

4'\x) - = v«)p')(c,).(x " x(o)). 

In order to extend Theorem [T| to the case of w 
we can write 



(11) 



w{x)-w{xm^w'{c).{x-xm, 

where w' is the Jacobian of the function w{x) and c is a 
point on the straight line between x and x(0) which may be 
different for different rows of w'{c) [10]. We may estimate 
the norm of w{x) — w;(x(0)) as follows: 



\\wix) - wixmm = ||u;'(c).(x-x(0))|| 

<lk'(c)||||(x-x(0))||. 



(12) 



The Lipschitz constant w may be estimated by maa;c||u''(c)|| 
where, || • || represents the Euclidean norm. 

We can apply the result of Theorem [T] to the nonlinear 
uncertain part of (fTOl l. Let us define a hyper rectangle 



X 

V 



Xi < Xi < Xj 

Vi <Vi < Vi 



(13) 



where Xiy ^^d Uj^ denote the lower bounds and Xi> ™d 
Vi denote the upper bounds on the new states and inputs 
respectively. For this purpose, the Jacobian of Wj-"'''(-) is 
found by differentiating it with respect to x and v at an 
arbitrary operating point — [x^ Pit)] for i — 1,2, ■ ■ ■ ,m 
and j = 1, 2, • ■ • , Ti where, x^ v ^ and p{t) G 51. Since 
we assume wp^(0, 0,p(t)) = 0, x(0) = 0, and v{0) = 0; 
wl''\-) can be written as follows: 



w'i'\x,v,Apit))^w'Y\c,,)-[x V Apit)] 
And then AW^(-) can be written as 



U), 



AW{-) = $ 



X 

Apit) 



(14) 



(15) 



where, 



$ = 



,l(l'2), 



Also, the bound on AW{-) can be obtained as follows: 

p = max||$||, i = l,2,---,m. (16) 

The bound in (fTST i is obtained by over bounding 

A. Linearized model with an Unstructured Uncertainty Rep- 
resentation 

In (fTsT l. C(.) is chosen such that it gives the maximum 
induced matrix norm on $. Once these bounds are obtained, 
we can write ( fTOl l in a suitable MIMO stochastic uncertain 
system form so that the minimax LQG control approach [7] 
can be utilized to design a tracking controller We define, 
Ci(^) = Ai(t)[(7ix(<) + Div(t)] G R", and W{t) = 



m w, 



r,TiT 



pm+n 



, where C, = Ci(0 + ^i' and wi is 



a disturbance input corresponding to Ap{t). Also W2 is a 
unity covariance noise input. We write linearized model as 
follows: 



X{t) = Ax{t) + Biv{t) + B2W{ty, 
z{t) ^ Cix{t) + Div{t); 
m = C2x{t)+D2W{t); 



(17) 



where. 



••■ 





P 
• • P 



Bi=B,B2=[ BiEi 0], Ci 
Di = ph, D2=[0 I ], ^[Ap{t)Ap{tf], 



Ci = E^^Ci, Di = E:[^Di and ||A|| < 1. Note that A{t) 
is m X m which satisfies the following stochastic uncertainty 
constraint condition. 



where 



E I" llCif <E r 
Jo Jo 

indicates the Euclidean norm.. 



(18) 



IV. MiNiMAX LQG Design 

The model developed in above section uses an uncertain 
diffeomorphism T{x,p{t)) which is unknown and hence 
any control system design using this model must contains 
a robust filter which able to estimates the uncertain states. 
Therefore, in this section we propose a minimax LQG design 
approach which uses a robust Kalman filter to estimates the 
elements of T{x, p{t)) and guarantees the stability and robust 
performance of the closed loop system. Here we present a 
summary of the minimax LQG design procedure. Interested 
readers are referred to [7] for more details on results and 
related proofs. 

The minimax LQG control problem [7] involves finding 
a controller which minimizes the maximum value of the 
following cost function: 



J : 



lim — E 

T->oo 2T 



(19) 



where i? > and G > 0. The maximum value of the 
cost is taken over all uncertainties satisfying the uncertainty 
constraint ( fTSl l. If we define a variable 



the cost function ( fT9] l can be written as follows: 



J: 



lim ( — )E 

T— »oo Al 



mvdt. 



(20) 



(21) 



The minimax optimal controller problem can now be solved 
by solving a scaled risk-sensitive control problem [7] which 
corresponds to a scaled iJoo control problem; e. g. se e [11 1. In 
this control problem the system is described bv jl7ll and (120b 
and the controller is to be constructed such that the closed 
loop system is stable and the transfer function from W{i) to 
satisfies the iJoc norm bound ||7Viii(j'w) || < IVw. The 
scaled risk-sensitive control problem considered here allows 
a tractable solution in terms of the following pair of iJoo 
type algebraic Riccati equations. 

- y^{Clv-^C-i - T~^R^)Y^ + BaC-f - Djr^^Ds)^^ = 0, (22) 

and 

- X^{BiG-^b'[ - T-^B2B^)X^ + (Br - T,G;'T^) = 0, (23) 

where, 

The solutions to both of the algebraic Riccati equations 
are required to satisfy the conditions Yoo > 0, X^o > 0, 
/ - r-iFoo^oo > and Rr - r^G-^r^ > 0. In order 
to solve the minimax LQG control problem, the parameter 
r > is chosen to minimize the cost bound (Wr) defined 
by 



Wr = tr[{TYCl + B.2Dl){D.2Dl)-^ 
X (rCaV + D2Bl)X{I - YXy^ + rVRr] 



(24) 



The minimax LQG controller H{s) has the following form: 



V = Kx, 



where, x £ R-" is the state of the controller and 

K = ^G-\B'[x^ + r'i:y, 

A, = A + BiK - B,C2 + r-\B2 - BcD2)B^X^. 
V. Air-breathing hypersonic flight vehicle 

EXAMPLE 

A. Vehicle Model 

In this section, we consider the same example as consid- 
ered in our previous work [5], [6]. The nonlinear model for 
the longitudinal dynamics of an AHFV is presented in [12]: 

T cos a — D . . L + T sin a q cos 7 

V= 5sm7, 7= -— , 

m mV V 

/l = Vsin7, d = Q - 7, Q = Myy/Iyy, (26) 



i = 1,2,3. 



Approximations to the forces and moments occurring in 
these equations are given as follows: 

Lf=:;g5'Ci(Q,5e,5c,Ari,Ar2), (27) 

D^q5'Ci3(Q,5e,5c,Ari,Ar2), (28) 

ztT + qScCuia, 5,, 5,, An, Ara), (29) 



My 



T ^ q[(PCTAa, An, Moo) + Cria, An, M^o, Ad)], (30) 

iV» « gCjvjQ,(5e,5c,An,Ar2], i = 1,2,3. (31) 

The coefficients obtained by fitting curves corresponding to 
these quantities are given as follows; here, we remove the 
function arguments for the sake of brevity: 

= eg a + C^" 5e + cl" 5c + C^"i An + C^"^ Ara + 0° , 

Cm = C%ia + cl} Se + C^J^c + C^/^ An + C^p Ara + Cli , 

Co = C^J^+^^i'' [a + An)' + c($'+'^"i'(a + An) + c'Js^, 



Ct 



+ C^;i An' + C^:;An+C«, 

Cn, = " + C^N- + <^N- + C^"' + C^"' + , 

(32) 

where n = [ni n2 n^]'^, and Ej G R^^'^ are vectors which 
describe the linear relationship Atj — Ejn for j — 1,2 [12]. 
The terms Moo and q are defined as follows: 



Mo- 



OS) 



(25) 



The nonlinear equations of motion in ( l26b have 11 states (in 
the vector x) in which there are five rigid body states; i.e., 
velocity V, altitude h, angle of attack a, flight path angle 7, 
and pitch rate Q and there are three vibrational modes which 
are represented by generalized modal coordinates Ui, where 
Ni is a generalized force. There are four inputs (in the vector 
u) and they are the diffuser-area-ratio Ad, the throttle setting 
or fuel equivalence ratio 0, the elevator deflection {6^), and 
the canard deflection Sc. For tracking control purposes we 
simpUfy the model is such a way that the scheme developed 



in Section |lll] can be used and the simplified model closely 
approximates the real model (see also, [13], [14]). Note that 
the effect of structural flexibility is entering into the system 
(|26] | through the forebody turn angle and aftbody vertex 
angle, Ari and At2 of the vehicle respectively. In the process 
of simplification, firstly, we remove all the flexible states 
fij for j — 1,2 from the CFM and consider the effect of 
flexibility in the model by considering Ati{x) and At2(x) 
as uncertain parameters. We simplify the forces and moment 
coefficients as follows: 

+ ACT,<t>{x, u), 
CL^Cia + d+ACiiu), 

Cm = Ct,a + [cIj - C%{--^)]6. + Clt + ACm{x, u), 

Cd = d^+^^'^\af + d^+^^'\a) + + ACd(x,M), 
Ct = C^" Aa + CTa + ~ M"" + C?. + ACt {x,u), (34) 

where ACa{-), ACi{-), ACt{-), ACtA') and ACm{-) rep- 
resent the uncertainties in their corresponding functions. 
Furthermore, in order to obtain full relative degree for the 
purpose of feedback linearization, we dynamically extend 
the system by introducing second order actuator dynamics 
(adding two more states cj) and 0) into the fuel equivalence 
ratio input as follows: 

(j) = ~2(bj„<j> - ui^(l> + bj^(l>c. (35) 

After this extension, the sum of the elements of vector 
relative degree will be equal to the order of the system n; 
i.e. n = 7 and thus satisfying one of the conditions for exact 
feedback linearization [8]. 

B. Robust Feedback Linearization of the Simplified Model 

The model obtained through the above simplification is 
still difficult to feedback linearize due to the presence of 
uncertainties in the system. We approach this problem by 
using the technique developed in Section |III] The outputs to 
be regulated are selected as the velocity V and the altitude h 
using two inputs, elevator deflection 5e and fuel equivalence 
ratio (f)c- Since 5c is a function of 5^; i.e related to 5e via an 
interconnect gain, we do not consider it as a separate input. 
Furthermore, we fix the diffuser area ratio Ad to be unity. 
This manipulation results in a 2-input and 2-output square 
system. The new simplified model consists of seven rigid 
states and two additional integral states as follows: 

x^[Vi V hi h a (I) ^ Q]^ , (36) 
where, 

Vi= I {V{t) - K)dr, hi= f {h{T) - hc)dT, 

and Vc and he are the desired command values for the 
velocity and altitude respectively. The uncertain parameter 
vector p e includes the vehicle inertial parameters, 
coupUng terms and the coefficients which appear in the 



force and moment approximations described previously and 
is given as follows: 

P^[C2 Cll cffACi ACd ACt 

ACm ACt^^I^gK^ (37) 

The model of the AHFV can be written in the form as 
follows: 

2 

= f{x,po) +^5fc(x,po)wfc + Af{x,p) 

k=l 

2 

+ '^^9k{x,p)uk; 

k=l 

y,{t) ^i^,{x,p), i = l,2 (38) 

where, Af{x,p) and Agk{x,p) are the uncertainty terms 
appearing in the corresponding functions. The control vector 
u and output vector y are defined as 

rrt 'J^ r^^ 'J~' 

u=[ui,u2] =[Se,(t>c] ,y=[yi,2/2] =[V,h] . 

We assume that p{t) G fi, where 57 is a compact convex set 
that represents the admissible range of variation of p{t) such 
that po lies in its interior In this study, a maximum variation 
of 10% of the nominal values has been considered. Thus, 
n = {p{t) e R9 I |0.9p,o| < \p^{t)\ < \l.lp,o\ for i = 
1, • • • ,9). It is worth mentioning the fact that there are no 
uncertainty terms exists in V, and h, we can write linearized 
input-output map for the original model ( |26] | using (|6]l as 
follows: 



- V " 












V 









AV- 


V 




Vl 




A V {x, u, p) 


h 







+ 





h 









Ah 


h 









A h (x, u, p) 






V2 




Ah'^{x,u,p) 



Corresponding uncertain diffeomorphisms for each system 
as in ^ which maps the new vectors ^ and rj respectively 
to the original vector x can be written as follows: 

^ = Ti{x,p{t),Vc), 7^ = T2{x,p{t),h,), (40) 

where, 

n{x,p{t),v.)=[ j^{v{T)-Vc)dT v-v, V vY, 

T2{x,p{t),K)='y jl{h(T)-h,)dT h~hc h h h], 

and Vc and he are the desired command values for the 
velocity and altitude respectively. Also, 

X = T{x,p(t),Vc,hc), (41) 

where X = [ Ci ?2 ^4 rji rj2 rjs rj4 rjs 1^, 

and T{x,p{t),Vc,hc) = [ Ti(x^p{t),Vc) T2{x,p{t),hc) f . 
And finally we can rewrite ( [39] l using the method given in 



Section Hill as follows: 







Vl 





f2 







Au)i(x, v,p{t)) 







(42) 



(43) 







Awi{x,v,p{t)) 





Furthermore, we can rewrite ( l42l i in the general form ( fTTI ) 
where, 



A: 






1 




















" 




" 











1 





































1 






















































1 




















1 











, B = 


























1 





































1 





































1 










































1 











10 X p 















































10 X p 







10 X p 






10 X p 



B2 = [0.1 X B Og], 



II A (i) 1 1 < 1 and x{t) £ is the state vector, v{t) = 
[vi V2]'^ E is the new control input vector and ({t) = 
A[Cix(i) + Div{t)] e is the uncertainty output. It is 
worth noting that the states ^1 , ^2, £,3, ?73, Vii ^75 can either 
be measured or constructed with the available hardware on 
board the aircraft. 

Finally, we design a minimax LQG controller as described 
in Section|IV]for AHFV velocity and altitude tracking control 
problem. The simulation result of the proposed controller 
with the original CFM (without simplification) in the pres- 
ence of time varying uncertainty is shown in Fig. [T] In the 
simulation certain velocity and altitude trajectories have been 
chosen to demonstrate the effectiveness of the approach. It is 
observed that the velocity and altitude tracking is achieved 
whUe control input and other states remain bounded. 

VI. Conclusion 

A robust nonlinear control scheme for an uncertain nonlin- 
ear system with time varying uncertainty is presented using 
robust feedback linearization and minimax linear quadratic 
Gaussian (LQG) methods. In the proposed method, a lin- 
earized uncertainty model is derived for the corresponding 
uncertain nonlinear system which is followed by a minimax 
LQG controller. The applicability of the scheme to a real 
world application is demonstrated by designing a robust 



tracking controller for an air-breathing hypersonic flight 
vehicle model with input coupling and flexible effects. The 
approach involves the linearization of a simplified curve 
fitted model using a robust feedback linearization method 
as the first step. In the second step, a velocity and altitude 
tracking controller is synthesized using the minimax LQG 
control design method. Simulation results with a large flight 
envelope simulation is also presented to demonstrate the ef- 
fectiveness of the scheme. The results show that the proposed 
method works very well under parameter uncertainties and 
give satisfactory results. Further, investigation of the results 
reveals that the minimax LQG based controller works well 
with parameter variations of up to 10% of their nominal 
values for which it is designed. 
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Fig. 1: Velocity and altitude reference tracking responses using the mean value approach. 



